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Abstract.  Single-zone  synchrotron  self-Compton  and  external  Compton  models  are  widely  used  to  explain  broad¬ 
band  Spectral  Energy  Distributions  (SEDs)  of  blazars  from  infrared  to  gamma-rays.  These  models  bear  obvious 
similarities  to  the  homogeneous  synchrotron  cloud  model  which  is  often  applied  to  explain  radio  emission  from 
individual  components  of  parsec-scale  radio  jets.  The  parsec-scale  core,  typically  the  brightest  and  most  compact 
feature  of  blazar  radio  jet,  could  be  the  source  of  high-energy  emission.  We  report  on  ongoing  work  to  test  this 
hypothesis  by  deriving  the  physical  properties  of  parsec-scale  radio  emitting  regions  of  twenty  bright  Fermi  blazars 
using  dedicated  5-43  GHz  VLBA  observations  and  comparing  these  parameters  to  results  of  SED  modeling. 


1.  Introduction 

Blazar  jets  are  known  to  radiate  across  the  entire  elec¬ 
tromagnetic  spectrum  (e.g.  IMarscher  20061 12010f).  Their 
radio  to  UV  (sometimes  -  X-ray)  emission  is  believed 
to  be  dominated  by  synchrotron  radiation  of  relativis¬ 
tic  electrons  while  radiation  at  higher  energies  could  be 
due  to  the  inverse  Compton  scattering  of  synchrotron 
photons  emitted  by  the  electrons  themselves  (the  syn¬ 
chrotron  self-Compton  process,  SSC,  IJones  et  al.  19741 
IGhisellini  fc  Maraschi  1989]l  and  photons  from  exter¬ 
nal  sources  (External  Compton,  EC,  ISikora  et  al.  19941 
IDermer  fc  Schlickeiscr  20021).  The  sources  of  the  external 
seed  photons  for  the  EC  process  include  the  accretion  disc, 
broad  line  region  (BLR)  clouds,  warm  dust  and  the  cosmic 
microwave  background  (CMB),  with  their  relative  contri¬ 
bution  vary  for  different  blazars. 

Despite  the  wide  acceptance  of  the  picture  outlined 
above,  the  exact  location,  geometry  and  physical  proper¬ 
ties  of  the  regions  responsible  for  blazar  emission  at  dif¬ 
ferent  bands  remain  largely  unknown.  Considerable  suc¬ 
cess  was  achieved  by  single-zone  models  (assuming  a  sin¬ 
gle  spherical  relativistically  moving  blob  of  magnetized 
plasma  which  emits  synchrotron  radiation  and  interacts 
with  an  ambient  photon  held)  in  explaining  blazar  emis¬ 
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sion  from  IR  to  y-ray  energies  (among  the  many  recent  ex¬ 
amples:  lAbdo  et  al.  20101  IFinke  fc  Dermer  20101  and  the 
references  in  Table  [1]).  A  signihcant  fraction  of  blazar 
radio  emission  is  known  to  originate  at  the  extended 
parsec-scale  jet  and  cannot  be  accounted  for  in  the  frame¬ 
work  of  a  single-zone  model.  However,  emission  from  in¬ 
dividual  jet  components  resolved  with  the  Very  Long 
Baseline  Interferometry  (VLBI)  can  usually  be  well  de¬ 
scribed  by  a  single  uniform  synchrotron  source  model 
dPacholczyk  1970]).  It  is  tempting  to  identify  one  of  the 
jet  components  resolved  by  VLBI  with  the  zone  responsi¬ 
ble  for  the  emission  at  higher  frequencies. 

In  this  contribution  we  describe  a  program  of  coordi¬ 
nated  multi-frequency  (4.6-43.2  GHz)  Very  Long  Baseline 
Array  (VLBA),  Swift  and  Fermi  observations  of  selected 
y-ray  bright  blazars  with  the  aim  to  identify  the  possible 
source  of  high-energy  radiation  within  the  parsec-scale  jet 
and  put  constraints  on  physical  properties  based  on  radio 
spectra. 

2.  Observations  and  analysis 

We  have  used  the  VLBA  ([Napier  1994D  to  observe  twenty 
blazars  that  were  expected  to  be  bright  y-ray  emit¬ 
ters,  prior  to  the  launch  of  Fermi-GST.  The  observations 
were  conducted  simultaneously  at  seven  frequencies  (4.6- 
43.2  GHz)  utilizing  the  array’s  unique  capability  of  fast 
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Table  1:  Comparison  of  the  emission  region  parameters 
estimated  from  the  VLBA  data  and  SED  modeling 


Multi-frequency  VLBA 

1  Published  SED  model 

results  (this  work) 

parameters 

Name 

RaZ  GHz  PvLBA 

[10^^  cm] 

Db 

^VLBA 

[G] 

-Dvar" 

Dsed 

Rsed 
[10^^  cm] 

PSED 

Bsed 

[G] 

Ref. 

AO  0235-1-16 

<  5800 

0.8 

<  11 

24.0 

B0528+134 

<  7600 

1.4 

<  1.2 

31.2 

S5  07163-714 

<  2100 

q:  = 

0.4'* 

10.9 

14 

40 

2.0 

1 

(1) 

OJ  248 

<  6300 

1.6 

<  23 

13.1 

OJ  287 

<  2800 

Q.  = 

o 

p. 

17.0 

W  Com 

<  1700 

0.8 

<  118 

1.2 

20 

3 

2.55 

0.35 

(2) 

3C  273 

<  1800 

2.0 

<  0.2 

17.0 

9 

20 

2 

12 

(3) 

3C  279 

<  4600 

1.4 

<  14 

24.0 

21.5 

25 

2.0 

1.8 

(4) 

PKS  B1510-089 

<  3900 

Q.  = 

0.2‘* 

16.7 

37 

18 

1.9 

0.09 

(5) 

4C  38.41 

<  3500 

1.0 

<  1.8 

21.5 

Mrk  501 

<  470 

1.6 

<  67 

f 

20 

1.03 

2 

0.3 

(6) 

NRAO  530 

<  6200 

1.6 

<  20 

10.7 

B1959+650 

<  510 

1.4 

<  300 

f 

18 

7.3 

2 

0.25 

m 

B2155-304 

<  2400 

1.0 

<  745 

f 

32 

150 

1.3 

0.018 

(8) 

BL  Lac 

<  590 

1.0 

<  3 

7.3 

3C  454.3 

<  6500 

q:  = 

0.8‘* 

33.2 

26 

15 

1.1 

5.4 

(9) 

B2344+514 

<680 

1.2 

<  57" 

^  p  is  the  power  law  index  in  the  electron  energy  distribution  N{E)  =  NoE~'^.  For 
the  optically  thin  part  of  the  synchrotron  spectrum  p  =  l  —  2a  where  a  is  defined  as 
St,  ~  ^  The  estimates  correspond  to  the  parsec-scale  radio  core.  The  value  is  in 

the  observer’s  frame.  The  homogeneous  synchrotron  source  model  is  not  applicable 
for  such  spectrum.  ®  The  variability  Doppler  factor  from|Hovatta  et  al.  (2009'y|used 
to  transform  the  observed  magnetic  field  strength  into  the  source  frame.  '  The  same 
Doppler  factor  as  in  the  corresponding  SED  model  was  adopted. 

References:  l:|Chen  et  al.  (2008)|  2;  [Acciari  et  al.  (2009)[  3:  |Pacciani  et  al.  (2009)| 
4:  [Giuliani  et  al.  (2009)[  5:  this  work;  6;  Anderhub  et  al.  (2009) [  7: 
[Tagliaferri  et  al.  (2008)|  ^  [Aharonian  et  al.  (2009)[  9:  Bonnoli  et  al.  (2010)[ 


frequency  switching  between  individual  VLBI  scans.  After 
the  initial  calibration  in  A  IPS  (IGreisen  1990|).  the  sources 
were  self-calibrated  and  imaged  independently  at  each 
frequency  using  Difmap  ([Shepherd  1994D.  This  software 
was  also  used  to  model  the  source  structure  in  the  uv- 
plane  and  derive  constraints  on  the  core  size  at  43.2  GHz 
(Table  H]). 

A  special  procedure  was  developed  to  improve  am¬ 
plitude  calibration  of  the  correlated  flux  density  result¬ 
ing  in  a  ^  5  %  accuracy  in  4.6-15.4  GHz  range  and  a 
^  10  %  accuracy  at  23.8  and  43.2  GHz.  Images  at  dif¬ 
ferent  frequencies  were  aligned  with  each  other  using  op¬ 
tically  thin  regions  of  the  jet  to  compensate  for  phase 
center  shifts  and  possible  frequency-dependent  core  po¬ 
sition  (“core  shift”,  ILobanov  19981  IKovalev  et  al.  20081). 
We  applied  a  spectrum  extraction  technique  described  by 


Sokolovsky  et  al.  (2010)  which  relies  on  the  positions  and 
flux  densities  of  individual  CLEAN  components  as  op¬ 
posed  to  the  use  of  restored  VLBI  images  in  order  to  sup¬ 
press  spectral  artifacts  arising  from  the  convolution  with 
a  Gaussian  beam. 

The  sources  were  observed  with  Swift  in  X-ray  and 
optical-UV  bands  typically  within  two  days  after  the 
VLBA  observation.  The  Swift  data  were  analyzed  us¬ 
ing  standard  tools  provided  in  the  HEAsoj^  package. 
Whenever  feasible,  standard  stars  within  the  UVOT  field 
of  view  were  utilized  to  improve  photometry  in  UBV 
bands.  The  Galactic  correction  was  applied  using  the  val- 


http : / /heasarc .nasa.gov/lheasoft/ 


Fig.  I:  Quasi-simultaneous  SED  of  PKS  BI5I0— 089 
constructed  using  observations  with  the  VLBA, 
Swift,  Fermi/ljAT,  NOT,  SMA  and  the  2.1  m  tele¬ 
scope  Guillermo  Haro.  The  violet  curve  (a)  represents  the 
combined  SED  model:  green  curve  (b)  is  the  synchrotron 
component,  red  (c)  is  the  SSC  component,  blue  (d)  is  the 
accretion  disk  and  brown  (e)  represents  EC  scattering  of 
the  disk  radiation  (see  also  Tabled]). 


ues  from  Schlegel  et  al.  (1998)  and  Kalberla  et  al.  (2005) 


for  the  optical-UV  and  X-ray  absorption  respectively. 

To  obtain  the  quasi-simultaneous  SED  of 
PKS  B1510— 089  (Fig.  [T])  discussed  below,  the  VLBA, 
Swift  and  Fermi/LAT  data  collected  during  two  days 
2009  April  09  -  10  were  complemented  by  UBVRci 
photometry  obtained  on  April  12  with  the  2.56  m  Nordic 
Optical  Telescope  (NOT)  and  infrared  JHKg  photometry 
obtained  on  April  07  and  17  with  the  2.1  m  telescope 
of  Guillermo  Haro  Observatory  in  Cananea,  Mexico. 
The  Submillimeter  Array  (SMA)  provided  flux  density 
measurements  at  1  mm  from  2009  April  05  and  14. 
These  observations  resulted  in  the  most  well-sampled 
quasi-simultaneous  SED  (containing  multi-frequency 
VLBI  data)  among  all  our  sources,  and  for  that  reason  it 
was  the  first  we  chose  to  construct  a  model  of. 

UBV  magnitudes  of  PKS  B1510— 089  observed  by  the 
Swift/WOT  are  systematically  (0U5-0U3)  brighter  then 
those  observed  by  the  NOT  two  days  later.  Since  both 
observations  were  calibrated  against  the  same  set  of  com¬ 
parison  stars  from  Gonzalez-Perez  et  al.  (2001)  the  differ¬ 
ence  can  be  interpreted  as  a  clear  sign  of  optical  inter-day 
variability. 


3.  Discussion  of  the  first  results 

Parsec-scale  radio  emission  of  all  the  observed  sources 
(Tabled])  is  dominated  by  a  bright  unresolved  core  (the  ap¬ 
parent  origin  of  the  jet,  see  a  discussion  bvIMarscher  2008|l 
which  exhibits  high-amplitude  flux  density  variability  (e.g. 
IVercellone  et  al.  2010|).  The  parsec-scale  core  is  a  natu¬ 
ral  candidate  to  be  directly  related  to  the  bright  variable 
emission  at  y-rays  ([Kovalev  et  al.  2009|1  and  other  bands. 
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1  10  100 
Frequency  (GHz) 


(a)  PKS  B1510-089 


(b)  OJ  248 


is  evidence  from  86  GHz  VLBI  data  that  H  ^  1  G  in  the 
mm-wavelength  core  of  3C  273  (|Savolainen  et  al.  20081). 

To  further  test  the  possible  relation  between  the  radio 
core  and  the  emitting  blob  implied  by  single-zone  SED 
models,  we  have  constructed  the  quasi-simultaneous  SED 
of  PKS  B1510— 089  presented  in  Eig.  [T]  The  use  of  the 
VLBI  core  flux  density  (instead  of  the  total  flux  density 
obtained  by  single-dish  radio  observations)  allows  us  to 
exclude  the  emission  coming  from  the  parsec-scale  jet  and 
larger  scales.  Our  ambition  was  to  check  how  this  will 
affect  the  SED  modeling. 


Eig.  2:  VLBA  spectra  of  core  regions  of  PKS  B1510— 089 
(a)  and  OJ  248  (b).  The  former  spectrum  is  flat  while  the 
latter  one  shows  a  prominent  synchrotron  self-absorption 
peak.  The  green  line  is  a  power  law  fit.  The  blue  curve  is 
a  homogeneous  synchrotron  source  model. 


The  multi-frequency  VLBA  observations  have  detected 
the  synchrotron  self-absorption  peak  in  core  regions  of 
most  of  the  observed  sources  with  the  exception  of 
B0716-H714,  OJ  287,  PKS  B1510-089  and  3C  454.3.  The 
latter  sources  show  inverted  spectra  in  the  4.6-43.2  GHz 
range,  characteristic  of  a  partially  optically-thick  inhomo¬ 
geneous  source.  These  sources  may  have  a  detectable  self¬ 
absorption  turnover  at  frequencies  >  43  GHz.  Pig.[2]com- 
pares  core  radio  spectra  of  PKS  B1510— 089  and  OJ  248. 

Eor  sources  where  the  self-absorption  peak  was  de¬ 
tected,  it  was  possible  to  put  upper  limits  on  the  mag¬ 
netic  field  strength  using  the  procedure  described  by 


Sokolovsky  et  al.  (2010)  The  source  intrinsic  values  were 


computed  using  variability  Doppler  factor  estimations 
from  Hovatta  et  al.  (2009)  The  resulting  liberal  upper 
limits  on  the  magnetic  field  strength  owe  to  poor  reso¬ 
lution  at  lower  observing  frequency  (4.6  GHz)  which  ulti¬ 
mately  limits  our  ability  to  unambiguously  extract  spec¬ 
trum  of  a  given  spatial  region. 

In  the  case  of  a  bright  isolated  component,  which  dom¬ 
inates  the  emission  across  the  whole  observed  frequency 
range,  it  is  possible  to  estimate  its  size  using  the  high¬ 
est  observing  frequency  (43.2  GHz)  instead  of  the  lowest 
one  and,  therefore,  obtain  much  tighter  constrains  on  the 
component  size  and  the  magnetic  field  strength.  The  ap¬ 
plicability  of  this  “isolated  component”  scenario  to  each 
individual  source  is  currently  being  investigated. 

Table  [T]  presents  constrains  on  size  (i?43  ghz),  mag¬ 
netic  field  strength  (Hvlba)  and  electron  energy  spec¬ 
tral  slope  (pvlba)  in  the  radio  core  region 


see 


Sokolovsky  et  al.  (2010)  for  a  more  detailed  discussion. 


These  values  are  compared  to  the  corresponding  parame¬ 
ters  of  SED  models  (i?sED,  Hsed,  Psed)  reported  in  the 
literature.  Doppler  factors  from  SED  models  (Dsed)  and 
those  assumed  for  Hvlba  computation  (Dvar)  are  also 
listed.  The  SED  model  parameters  are  well  within  the  con¬ 
straints  set  by  the  multi-frequency  VLBA  data,  except  for 
3C  273  which  shows  a  discrepancy  in  Hvlba-JJsed  that 
may  be  due  to  insufficient  resolution  of  our  data  set.  There 


Table  2:  PKS  B1510— 089  SED  model  parameters 


Minimum  e~  Lorentz  factor 
Brake  Lorentz  factor 
Maximum  e~  Lorentz  factor 
e~  energy  slope  below  yurk 
e~  energy  slope  above  yurk 
Doppler  factor 
Bulk  Lorentz  factor 
Magnetic  field  strength 
Variability  time 
Blob  radius 

Jet  power  (magnetic  field) 
Jet  power  (electrons) 

Black  hole  mass 
Accretion  efficiency 
Eddington  ratio 
Blob  distance  from  disk 


ymin  =  30 

ybA  =  1.9  X  10® 

ymax  =  1.0  X  10® 

Pi  =  1.9 
P2  =  3.9 
D  =  37 
r  =  37 
5  =  0.09  G 
tvar  =  2.16  X  lO"'  sec 
5=1.8  X  10®®  cm 
Pjet.  B  =  2.8  X  10®®  erg/s 
■Pjet  e-  =  7.5  X  10®®  erg/s 
IUbh  =  1  X  IO^ALq 
y=l/12 
Ldisk/LEdd  =  0.06 

r  =  5.1  X  lO®®”  cm 


Contribution  from  a  hot  thermal  component,  proba¬ 
bly  the  accretion  disk,  is  evident  in  the  optical  and  UV 
bands.  The  SED  in  the  IR  to  y-ray  range  can  be  ex¬ 
plained  by  the  single  zone  EC  model  with  the  Compton 
scattered  accretion  disk  photons  being  responsible  for  the 
bulk  of  the  observed  y-ray  emission.  In  this  model,  the 
blob  is  close  enough  to  the  accretion  disk  so  the  accre¬ 
tion  disk  photons  can  reach  the  blob  directly,  with  no 
need  to  be  scattered  back  from  BLR  clouds  or  dust.  Note, 
that  due  to  relativistic  aberration  some  of  the  disk  pho¬ 
tons  will  be  coming  to  the  blob  (nearly)  head-on.  The 
photons  will  primarily  come  from  a  radius  of  the  disk 
of  \^r  where  r  is  the  blob’s  distance  from  the  black 
hole  (jPermer  fc  Schlickeiser  2002j) .  The  blob  size  was  con¬ 
strained  by  the  observed  y-ray  variability  time  scale  of  the 
order  of  a  few  hours  (|Abdo  et  al.  2010)) .  Parameters  of  the 
model  are  summarized  in  Table  [H  We  note  that  the  model 
fit  is  not  unique,  single-zone  models  with  a  different  source 
of  seed  photons  could  probably  also  produce  an  acceptable 
fit  (cf.  lAbdo  et  al.  2dT0)). 

The  constructed  model  (Fig.  [T|)  represents  well  the  IR 
to  y-ray  data,  however  it  dramatically  underpredicts  the 
radio  emission  because  of  the  synchrotron  self-absorption 
occurring  in  the  small  blob  at  a  high  frequency.  A  sig¬ 
nificantly  larger  blob  size  would  be  inconsistent  with 
the  observed  variability  timescale.  Even  if  the  variabil¬ 
ity  timescale  argument  were  dismissed,  a  single-zone  SED 
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model  with  a  larger  blob  would  have  difficulties  in  ex¬ 
plaining  the  observed  hard  radio  spectrum  which  is  incon¬ 
sistent  with  the  softer  X-ray  spectrum  or  optically  thick 
synchrotron  radiation  from  a  uniform  source.  One  way  to 
overcome  this  difficulty  would  be  to  introduce  a  second 
brake  in  the  electron  energy  spectrum. 

Alternatively,  one  could  abandon  the  attempt  to  ex¬ 
plain  the  whole  SED  with  a  single-zone  model  and  assume 
that  even  the  radio  emission  observed  from  the  cm-mm 
band  core  is  coming  from  a  larger  (Ased  <  R  <  R43  ghz) 
structure  downstream  of  the  blob  responsible  for  the  IR- 
to-7-ray  emission.  This  structure  (which  cannot  be  ob¬ 
served  separately  from  the  blob  due  to  limited  resolu¬ 
tion  of  the  available  VLBI  data)  could  be  a  smooth  jet 
dBlandford  &:  Konigl  1979[)  or  a  number  of  distinct  jet 
components  ((Marscher  198011  -  perhaps  the  blobs  which 
contributed  to  high-energy  emission  earlier...  One  could 
argue  that  currently  this  extended  region  does  not  con¬ 
tribute  significantly  to  the  optical  and  7-ray  emission  on 
the  basis  of  the  observed  short  timescale  variability  in 
these  bands.  A  more  detailed  SED  modeling  is  needed  to 
test  this  scenario. 

4.  Summary 

The  program  of  coordinated  multi-frequency  (4.6- 
43.2  GHz)  VLBA,  Swift  a.nd  Eermi  observations  of  selected 
7-ray  bright  blazars  is  described  and  the  first  analysis  re¬ 
sults  are  discussed.  Constraints  on  the  size  (it!43  GHz)  and 
magnetic  field  strength  (Bvlba)  in  the  parsec-scale  core 
region  derived  from  the  VLBA  observations  are  consistent 
with  the  values  (i?sED,  Bsed)  usually  assumed  in  single¬ 
zone  SED  models  (Table[T|).  However,  the  single-zone  SED 
model  which  we  tested  for  PKS  B1510— 089  has  difficul¬ 
ties  in  explaining  the  observed  level  of  core  radio  emission 
while  being  consistent  with  the  observed  7-ray  and  opti¬ 
cal  variability  timescales.  The  multi-frequency  VLBA  and 
SED  data  analysis  for  the  whole  sample  is  in  progress.  The 
SED  data  on  PKS  B1510— 089  discussed  here  are  available 
from  KVS  upon  request. 
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